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Abstract: In pursuing the intriguing resemblance of the Einstein equations to 
thermodynamic equations, most sharply seen in systems possessing horizons, we 
suggest that eternal inflation of the stochastic type may be a fruitful phenomenon 
to explore. We develop a thermodynamic first law for quasi-de Sitter space, valid 
on the horizon of a single observer's Hubble patch and explore consistancy with 
previous proposals for horizons of various types in dynamic and static situations. 
We use this framework to demonstrate that for the local observer fluctuations of 
the type necessary for stochastic eternal inflation fall within the regime where the 
thermodynamic approach is believed to apply. This scenario is interesting because 
of suggestive parallels with black hole evaporation. 

Keywords: Effective Field Theory, Cosmology, Inflation, Thermodynamics. 



Contents 



1. Introduction 



1 



2. Defining the horizon 



4 



3. First law of thermodynamics for the apparent horizon 

3.1 Energy flow and surface gravity 

3.2 Defining 5E 

3.3 Defining surface gravity 

3.4 Einstein equations as the first law 



8 

8 
9 
10 
12 



4. Including Perturbations 

4.1 Perturbing the spherically symmetric metric 

4.2 Comparing the 'thermodynamic' equation 



14 

14 
16 



5. Eternal Inflation 



16 



6. Discussion 



18 



A. Surface gravity 



19 



A.l Review of surface gravity for a static spherically symmetric space-time 19 
A. 2 Surface gravity in time dependent spherically symmetric spacetimes 22 



1. Introduction 

Black hole thermodynamics [1, 2], from the most utilitarian perspective, provides 
a practical shortcut for describing the black hole evaporation process. In order to 
predict what happens to an isolated black hole, one should in principle work out the 
renormalized energy-momentum tensor of the known quantum fields on a black hole 
solution and plug it back into the Einstein equations. This is known to be a very 
complicated calculation, involving non-local expressions and, in four dimensions, 
semiclassical equations of fourth order. Two-dimensional examples of the precise 
calculation are well studied, largely following the work of Callan, Giddings, Harvey 
and Strominger [3], but the four dimensional case is so far much less tractable. 
However, several interesting features of black hole evaporation can be understood 
from the thermodynamic formulation alone, and in this respect it offers a great 
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simplification [4]. The rate of emission of a black hole is that of a black body of 
temperature T = 1/8ttM with grey body factors obtainable in perturbation theory 
for each species. We can thus sketch the geometry of an evaporating black hole as an 
approximate Schwarzschild solution with adiabatically varying mass M ~ — 1/M 2 
and easily estimate the evaporation time t cv ~ G 2 M 3 ~ RS, where G is Newton's 
constant, R is the radius of the black hole and S is its entropy. 

What are the limits of applicability of the above description? Naively, the semi- 
classical treatment is accurate enough for macroscopic black holes whose mass is 
large in units of the Planck mass, M 3> Mpi, and on scales longer than the Planck 
scale. That is, the scenario seems likely to be under control as long as the back- 
reaction of the quantized fields is small compared to the typical curvature at the 
horizon. However, several arguments, nicely summarized in [5], indicate that the 
naive perturbative low-energy approach is only valid on scales parametrically longer 
than the Planck scale and set by t cv ; this is a statement of the black hole information 
puzzle. 

It is tempting to speculate that the thermodynamic story can be extended to 
systems beyond black holes, where it might similarly capture aspects of quantum 
gravity independent of the microphysics at work on very small scales. Some au- 
thors have gone even further and suggested the ambitious program of deriving the 
Einstein equations from thermodynamics. Jacobson extended the first law of black 
hole thermodynamics to any local Rindler horizon and showed that dQ = TdS, with 
suitable definitions of those quantities, contains the same information as the Einstein 
equation [6] (see also [7, 8, 9]). Reference [10] contrasts some of these approaches. 
Hayward [11] and more recently Cai and collaborators have considered spherically 
symmetric spaces, including FRW [12, 13, 14]. Multiple fluids in a cosmological 
setting have been considered in [15]. Frolov and Kofman [16] considered a quasi-de 
Sitter space and showed that the evolution of the (perturbed) apparent horizon can 
also be reformulated as the first law of thermodynamics. Inspired by more recent 
speculations by Verlinde [17] and subsequent investigations (e.g. [18]) it has been 
shown [19] that a thermodynamic relation can also be applied on a time-like screen, 
once a few geometrical constraints are met. At the current level, most of the attempts 
to derive the Einstein equations in fact invoke some knowledge from them (for ex- 
ample in starting with a metric that is a vacuum solution or in choosing a definition 
of mass or energy flux justified by relating the stress-energy tensor to geometry) and 
so it is difficult to judge the success of that approach so far. 

From a more conservative perspective, the first step in exploring the thermo- 
dynamic framework is to look for cases where geometric quantities take properties 
similar to those we associate with temperature and entropy. Since the evaporation is 
the key feature of interest for the black hole case, it may be particularly interesting to 
consider thermodynamic descriptions of other systems where the back reaction of the 
quantized fields plays a relevant role. Beyond the computationally practical reason 
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highlighted in the first paragraph, such systems may provide a means to search for 
and clarify puzzles similar to the black hole information story. Eternal inflation of 
the stochastic type is a well-known example of such a system and so it is worthwhile 
to consider whether a thermodynamic treatment can be applied there. 

Our work will be similar to the early black hole work, rather than Jacobson's 
more general approach: the analysis uses a specific class of solutions (spherically 
symmetric, assuming adiabatic evolution). We will not try to derive the Einstein 
equations but rather understand whether a thermodynamic interpretation can apply 
in the stochastic eternal regime. 

For a slow-roll inflation scenario we can exploit the nearness to de Sitter and 
assume that the temperature associated with the Hubble horizon is T ~ H , where 
Hq is the Hubble parameter at the time of interest. For a spectator scalar field in an 
exact de Sitter background, one can calculate the contribution from the fluctuations 
of the field to the background renormalized energy momentum tensor, e.g. (T^) rcn ~ 
Hq ~ T 4 [20]. This is always a negligible contribution to the total energy density 
Ptot ~ HqMpy. In scenarios close to de Sitter considering a single Hubble patch 
this calculation is expected to hold qualitatively but it is not the only effect of the 
fluctuations in the eternally inflating regime. There the change in field position due 
to the fluctuations is as large as the classical rolling of the field down the potential 
in a Hubble time. The fluctuations qualitatively change the evolution of the system 
even though the Hubble parameter itself may (locally) remain nearly constant. In 
practice one takes the classical position of the field on the potential to be adjusted by 
the combination of classical rolling and quantum fluctuations, where this adjustment 
is homogeneous on a Hubble scale and is performed in steps of order the Hubble 
time. 

Intuitively it seems reasonable to expect that the phenomena of stochastic eter- 
nal inflation, like black hole evaporation, does not depend on a precise formulation 
of quantum gravity. In order to better define the domain of validity of effective field 
theory in this context we might try to explicitly suppress the quantum gravity ef- 
fects (which in the black hole evaporation case are responsible for the recovery of 
information) and then ask about the validity of our description. The simplest de- 
coupling limit, M p — > oo, is too trivial as it makes gravity oblivious of the matter 
fields' dynamics and kills the scalar fluctuations that are responsible for the chaotic 
inflationary regime. But because there is an extra parameter in inflation that con- 
trols the departure from equilibrium (the slow-roll parameter e = —jp), we might 
hope to decouple UV gravitational effects while still maintaining a dynamical system. 
Naively, together with sending Mpi — > oo, we want to send the combination H 2 /e to 
infinity in such a way that the amplitude of curvature perturbations 
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remains constant (and such that we return neither to the flat space limit, nor the limit 
where the horizon is too small for quantum modes to created). Tensor modes and 
graviton loops will be suppressed, and we can conclude that any remaining puzzles 
are related to IR effects. Although we leave a rigorous analysis of this decoupling 
limit for further work, it seems very reasonable to expect that this limit or some 
variant is possible. If the stochastic eternal inflation scenario has a description at sub- 
Planckian scales it can be a useful test case for the thermodynamic story, analogous 
to the evaporating black hole. 

To that end, we will follow up on a relatively conservative extension of the 
thermodynamic ideas as considered by Frolov and Kofman [16], where we allow a 
quasi de Sitter space and small metric fluctuations, requiring spherical symmetry 
for simplicity. We will first show that one can modify their approach to include 
contributions from pressure to the first law and that it then yields a temperature 
consistent with Hawking temperature. We demonstrate that, regardless of what the 
deep implications may be, the thermodynamic relationship as proposed holds on 
the horizon of a single observer in near-equilibrium systems during both slow-roll 
and stochastic eternal inflation. For slow-roll inflation, Frolov and Kofman gave 
a treatment that included the first order fluctuations of the metric and inflaton 
field. However, they left the eternal inflating regime out on two grounds: first, that 
increasing Hubble parameter is incompatible with the classical Einstein equations; 
second, that the large variation in the Hubble parameter on super-horizon scales 
indicates that an adiabatic approach is not sufficient. However, the first objection 
is no obstacle in the black hole evaporation scenario, where the back-reaction of 
quantum fluctuations also does not behave as a classical matter source. The second 
objection is removed if we take a local perspective and follow a single Hubble patch: 
for a local observer the eternal inflation regime is (at least for a short time, and in 
the slow-roll scenario) perturbative. 

The plan of the paper is as follows. In the next section we review the definition of 
the apparent horizon. In Section 3 we set up the first law of thermodynamics in the 
quasi-de Sitter context and evaluate that equation for slow-roll inflation. In Section 
4 we include fluctuations. In Section 5 we demonstrate that the thermodynamic 
approach finds the same success and difficulty in the eternal inflating regime as 
the standard general relativistic approach. In the final section we discuss some 
implications of this result and conclude. 

2. Defining the horizon 

The standard ingredient for applying the thermodynamic relation dQ = TdS = 
dE + PdV to spacetime is the presence of an event horizon, which hides information 
and can thus be associated with an entropy S. In this context the horizon is inevitably 
global and observer-dependent; it requires the knowledge of the entire spacetime and 
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that of the observer's world-line in order to be defined. For practical purposes it is 
useful to look for a different type of horizon, and deal with quantities that can be 
locally characterized. A review of these issues can be found in [7]. 

In the literature different types of surfaces have been proposed that can be com- 
patible with the notion of entropy. In the vicinity of any spacetime point there exist 
congruences of null rays (with locally vanishing expansion) that form a local Rindler 
horizon for a class of accelerating observers. Jacobson [6] applies the first law of 
thermodynamics to such null rays. He considers the variation of area spanned by the 
congruence of null-rays and relates it with the flux of matter through the Rindler 
horizon. Since the congruence is initially non expanding, dS is given by integrating 
the logarithmic derivative of the expansion which is related to Ricci tensor through 
the Raychaudhuri equation. In addition, Jacobson's association of this surface with 
entropy is based on the notion of entanglement entropy. 

Frolov and Kofman [16], on the other hand, focus their attention on the apparent 
horizon in spherically symmetric spaces. Their approach is slightly different because 
they follow the apparent horizon (which is generally a non-null surface) and define 
dS through the variation of its area. The definition of the apparent horizon relies on 
the existence of a marginally trapped surface: a two-dimensional space-like surface 
orthogonal to a congruence of null rays with zero expansion. A marginally trapped 
surface that is also closed generally indicates the presence of strong gravitational 
effects. The set of marginally trapped surfaces forms an apparent horizon. We will 
follow Frolov and Kofman and use this surface to define a thermodynamic analogy 
applicable to inflationary spacetimes (this is also similar to what is proposed for 
dynamical black holes, e.g. [21, 22, 23, 24]). The association of entropy proportional 
to the area of this surface comes by analogy with the black hole (for de Sitter, the 
Gibbons-Hawking entropy [25]). Currently, the most compelling direction to explain 
all the area/entropy laws may be the notion of entanglement entropy, but no exact 
calculation has proven this. 

In an expanding cosmology, in order to find a closed marginally (anti-) trapped 
surface, one has to consider larger and larger spheres, until the expansion of the 
Universe counterbalances the tendency of the ingoing light rays to "converge" . More 
explicitly, consider a flat FRW metric 



Radial light rays satisfy the equation dr = ±dt/ a, where a plus sign holds for outgoing 
and a minus for ingoing rays respectively. A beam of light rays leaving a sphere of 
constant r describes two spheres ("+" and "— ") of areas A±(t) = 47ra 2 (t)r±(t), where 
r±(t) are the two solutions for the radial null rays. The derivative of the area along 
the two beams, 



ds 2 = -dt 2 + a 2 (t)(dr 2 + r 2 dtt). 



(2.1) 




(2.2) 
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Figure 1: We plot some relevant hypersurfaces of an inflationary spacetime. For illustra- 
tive purposes we take e = 0.4. The axes represent the comoving radius r and conformal 
time r; radial light rays are lines at 45 degrees. One of them is the future event horizon 
(EH) for the observer sitting at r = 0, defined as the past light cone of the observer at fu- 
ture infinity [26]. What happens outside the EH has no influence on the observer. We also 
plot the apparent horizon (AH), defined by r = (aH)" 1 . The dashed curves are examples 
of 2-dimensional spheres of constant area. When these surfaces intersect the apparent hori- 
zon their tangents are at 45 degrees with respect to the coordinate axes. In other words, 
the AH is the place where incoming light rays have (instantaneously) zero expansion. The 
vectors £ and £, on the AH, are parallel to the AH and to the null direction respectively. As 
summarized in the text, Jacobson's analysis [6] applies to a null surface, relating the flux 
of matter through it with the expansion of the null generators (i.e. along £ in this context). 
Here, as in [16, 12], we consider the system contained within the AH i.e. we study the 
matter flux through the AH and relate it with the expansion along the £ direction. The 
apparent horizon, event horizon, and the surface p = const. = H" 1 all coincide in the de 
Sitter limit e = 0. 



has a familiar behavior: at small radii (ra < H^ 1 ) the light rays going inward tend to 
contract, the ones going outward tend to expand. However, sufficiently large spheres 
"feel" the expansion of the Universe. For r > (aH)~ l both classes of light rays are 
expanding and the two-dimensional surface is called (anti-)trapped. At any time, the 
two-dimensional surface with r = (aH)~ l is thus marginally trapped and is referred 
to as apparent horizon. The union of all the marginally trapped surfaces, i.e. the 
three dimensional surface, E : r(t) = (a(t)H(t)) , is called a trapping horizon. 
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Note that a trapping horizon is generally not a null surface; for e > it is time- 
like. For evolving black holes, a great deal of work has gone in to defining local notions 
of horizons, especially by Ashtekar and collaborators, but this has focused on null 
and space-like surfaces [22, 24] which are relevant for stationary and accreting black 
holes. Interestingly, it is the third case (time-like surfaces) that is most needed for 
classically evolving quasi-de Sitter space. For more details about trapping horizons 
we refer readers to Appendix A. 2. As is common in literature, in what follows we 
will refer to both two dimensional apparent horizon and three dimensional trapping 
horizon, E, as the apparent horizon. 

Frolov and Kofman [16] associate the entropy S = A /AG to the area A, corre- 
sponding to two dimensional apparent horizon at anytime. For an inflationary FRW 
universe with a small and constant slow roll parameter e = — -J^, it is easy to show 
that the ratio of the area of the apparent horizon to that of the future event horizon 
is close to one (see Figure 1) 

^ = l-2e. (2.3) 

■"■eve 

Since we will eventually consider perturbations of the metric (2.1) it is useful to 
work with a more general metric ansatz. We will maintain spherical symmetry of the 
background in order to make plausibly well defined identifications between geometric 
and thermodynamic quantities. For this purpose consider the metric 

ds 2 = g ab dx a dx b + p 2 (t, r)dQ , (2.4) 

where a, b = 0, 1 label the two coordinates t and r, and p is a generic function of the 
latter. We will use greek indices a, f3 for four-dimensional quantities. In our notation 
V represents the covariant derivative in 1 + 1 dimensions (with respect to g a b) and 
we use the symbol (;) to indicate the covariant derivative in 3 + 1 dimensions. 

The variation of the area of a surface along a trajectory x a {\) is given by 
dA = k a \7 a Ad\, where k a = dx d [ X ^ and A represents displacement along the tra- 
jectory that can be taken to be the time coordinate for time-like or null trajectories. 
By generalizing the argument given in Eq.(2.2) we find the trapped surface for a 
spherically symmetric metric by considering 

dA dx a (t)„ . , 

* = 8 *'^r v «" = °- (2 - 5) 

where dx a is null. Since V a p is also a temporal-radial vector, this equation implies 
that V a p should be a null vector. The condition for a marginally trapped surface in 
the spherically symmetric case is then: 

(Vp) 2 = (2.6) 

Setting the position of apparent horizon by / = (Vp) 2 = 0, V a / is normal to 
apparent horizon, while 

C = e ab V b f (2.7) 
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may be used to represent the (non-angular) tangential vector along this surface. Here 
e ab is the two-dimensional (t, r) Levi-Civita tensor and ( a for indices other than (t, r) 
are zero. 

For cosmological scenarios ( a is time-like, while for exact de Sitter space, ( a 
becomes a null vector and so is parallel and normal to the usual time-like Killing 
vector at the de Sitter event horizon (which is also the apparent horizon). Now 
calculating k a along the non-angular generators of the surface, since both ( a (\) an d 
k a (X) are tangential to this trajectory they must be proportional for all values of A. 
For simplicity we assume ( a = B(X)k a and in the end -B(A) will be irrelevant for our 
analysis. Therefore, the variation of entropy is related to the variation of the area of 
the apparent horizon by: 



In analogy with the black hole case, this is the entropy associated with things outside 
the horizon measured by an observer inside. In the next section we will consider 
appropriate definitions of mass, pressure and temperature to write the first law. 

3. First law of thermodynamics for the apparent horizon 

Here we will obtain a first law applicable during inflation which can also incorporate 
perturbations for spherically symmetric solutions. This derivation differs slightly 
from what was proposed in [16] in a way which makes the quantities identified with 
temperature, mass, and pressure more comforting but makes no difference in the 
actual calculation: the inputs into both their calculation and ours are geometric 
identities and the Einstein equations. The intermediate step of labeling thermody- 
namic quantities may give us some physical insight, but at this level is nothing more 
than a labeling. Our definitions are similar to those used for spherically symmetric 
spaces (without fluctuations) by Hayward [11] and more recently Cai et. al. [14], 
although we are careful to distinguish signs and the physical interpretation for de 
Sitter compared to black holes. In addition we follow how the slow-roll parameter e 
controls the departure from equilibrium and the corrections to the static definitions 
of geometric and thermodynamic quantities. As long as e < 1 the system is near- 
equilibrium and deviations from exact de Sitter provide only higher order corrections 
to the thermodynamic equation we will find. However it may be important to keep 
track of the e corrections since the system of real interest here cannot have e = 
(in that case there is only a trivial limit for decoupling gravity and the scenario may 
lose much of its power as a test case). 

3.1 Energy flow and surface gravity 

To make the analogy between the Einstein equations and the first law, we need 
to introduce quantities that will play the role of temperature and energy. Precise, 
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globally well defined notions of mass and surface gravity (temperature) are only 
possible for static, asymptotically flat spaces, although a great deal of work has been 
done to define quantities that are useful for astrophysical, evolving black holes. Here 
we will exploit spherical symmetry in order to define local quantities, and we will 
ensure that the standard results are recovered for static Schwarzschild and de Sitter 
spacetimes. 

For a spherically symmetric metric, given in Eq.(2.4), the relation between the 
four dimensional Ricci scalar and that derived with the two dimensional metric 
g a b is 

^R = R-^ + -\l-(Vp) 2 }, (3.1) 
P P 

where Op is defined from the two-dimensional metric. By integration over the angular 
variable we obtain the Einstein Hilbert action for a spherically symmetric space, 



Seh = 4:71 J d xy— g 
and the (a b) Einstein equations 



16nG 



(p 2 R + 2(Vp) 2 + 2) + p 2 C 

matter 



(3.2) 



2p- 1 [g ab n P - V fe V Q p] + p- 2 gab[{Vp) 2 - 1] = 8nGT ab . (3.3) 
3.2 Defining 5E 

To associate an appropriate quantity to the change in energy, 5E, we make use of 
the fact that for spherically symmetric but not necessarily static systems there still 
exist conserved currents [21] based on which one can define the flux transfer across 
a hypersurface. In particular we will use a definition of 5E which is related through 
the Einstein equations to 5M where M is the Misner-Sharp mass, a well-defined local 
quantity for spherically symmetric metrics. 

While in time-dependent space-time there is no Killing vector for time transla- 
tion, Kodama has shown [27] (also see [28]) that for spherically symmetric metrics 
there is a divergence free vector £ a such that, 

C = -e ab V hP , (3.4) 

and £ a = for indices other than (t, r). It is easy to show this vector satisfies = 
and is null on the surface of the apparent horizon. In addition Kodama proved that 
G 01 ^^ also corresponds to a conserved current [27]. Applying the Einstein equations 
to this expression we can define a current: 

j a = ~T a % (3.5) 

which is also covariantly conserved, j a ; a = 0, and the minus sign ensures the flow 
for quasi-de Sitter is future directed. This enables us to calculate the energy flux 
through the apparent horizon, 

5E = ] a dZ a . (3.6) 
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dH a denotes the 3-volume element of the apparent horizon and it can be written in 
terms of ( a : 

dE a = 4np 2 k b e ab d\ 

= (B- l d\)4Trp 2 ( b e ab , (3.7) 

while dH a = for indices other than (t,r). 

Although we will consider the current above to be the definition of 5E, it is 
important to emphasize how the choice is motivated. The original current of Kodama 
was geometric, and the current j a above is related through the Einstein equations to 
the same well-defined quasi-local quantity for spherically symmetric metrics [28], 

f = ^ 2 e ba V b M, (3.8) 

where M is defined by 

M=^[l-(Vp) 2 ]. (3.9) 

In the literature M is referred to as Hawking-Israel/Hernandez/ Misner- Sharp mass 
and a thorough discussion of its properties can be found in [21]. Now combining 
Equations (3.6), (3.7), (3.8) and (2.7), we obtain 

5E = (B- 1 dA)C°V tt M (3.10) 

Taking a derivative of above mass and using the (two-dimensional) Einstein equa- 
tions, yields that 5E is also related to stress energy tensor in the following way: 

5E = ^(B-'dX) [(1 - Vp 2 )V aP - pV a (Vp) 2 ] C 

= {B- l d\)^p 2 {T ab - Tg ab )V b p C . (3.11) 

3.3 Defining surface gravity 

The other key ingredient we need is a prescription for the surface gravity, k. The 
standard geometric definition is only applicable in a stationary spacetime in the 
presence of a Killing horizon and one can propose various strategies to extend it to 
the dynamical apparent horizon (the marginally trapped surface). In what follows 
we take the approach proposed by Hayward [11] for spherical symmetric metrics 
which produces expected values for known stationary cases. The new definition can 
be expressed as, 

£°V [a &] = K6 (3.12) 

or equivalently as, 

«=i|Dp|. (3.13) 
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For details on why these expressions agree, even up to the sign, and also how they 
relate to standard definitions of surface gravity we refer the readers to Appendix 
A. It is difficult to be rigorous, so we content ourselves here with the fact that this 
definition gives k = H for static de Sitter, and k = for Schwarzschild, and that 
our results won't be sensitive to the distinction between this definition and others 
that reduce to the correct answer in the static case. As discussed in the appendix, 
there are several proposals for defining surface gravity for non-static cases. They 
all agree in the static limit but may differ in corrections proportional to e. Since 
our thermodynamic expression will already be first order in e, these differences in 
definition do not affect our basic result. 

The crucial element for considering event horizons as thermodynamic objects is 
Hawking radiation [29, 30, 31, 25, 32], which has a thermal spectrum with temper- 
ature T = j-. In exact de Sitter space the temperature of quantum fluctuations 
is related by the same formula to the surface gravity k^s = H that any comoving 
observers associates with their de Sitter horizon. Defining temperature for a dynam- 
ical horizon is a more challenging issue. In highly dynamical situations the system 
does not have time to reach an equilibrium and the canonical notion of temperature 
does not exist. This implies that even though one may still interpret which term 
represents the heat absorption in a meaningful way, there is no clear way to dis- 
tinguish between entropy and temperature. Similarly for highly dynamical horizons 
one does not expect the existence of a quantity that can be directly interpreted as 
an instantaneous physical temperature of some radiation. However, if we consider 
slowly evolving scenarios (small departures from de Sitter) it seems reasonable to 
associate a temperature to the local geometry at any particular time. Therefore, 
for quasi-de Sitter space with e = < 1, we expect that any corrections to the 
temperature be proportional to e, that is T = ¥-{1 + 0(e)). One may be able to 
make more precise statements about the adiabaticity conditions following the work 



Note that in the laws of equilibrium thermodynamics it is also usually considered 
that transitions happen from one nearly time independent state to another one. This 
implies that in infinitesimal physical processes the change in surface gravity in the 
term that represents TdS is negligible the same way that change in the temperature 
is ignored in standard thermodynamics. Inferring temperature as T = ^- using the 
definition of Equation (3.13) we obtain 



which agrees with our intuition. 

The temperature associated to the apparent horizon by Frolov and Kofman [16] 
was different as they interpreted a different geometric quantity as k, based on a 



of [33] or [34]. 




(3.14) 
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guessed analogy with the black hole geometry: 

n FK = ^-[l-{Vpf] . (3.15) 

They wrote a geometric equation that looked similar to the first law (with zero pres- 
sure) where this quantity, Kfk, appeared as a constant of proportionality between dE 
and dS. However, this choice evaluates to H/4tt on the horizon, half of the standard 
Hawking temperature. We will write a similar geometric equation but identify the 
temperature as the geometric quantity given in Eq.(3.13). The difference between our 
results and theirs underscores a primary difficulty in the thermodynamic interpre- 
tation: there are many mathematically correct combinations of geometric identities 
and the Einstein equations for inflating space times, but without the additional input 
(e.g. calculating the Hawking temperature) it is difficult to unambiguously determine 
which quantities (if any) have truly meaningful thermodynamic interpretations. 

3.4 Einstein equations as the first law 

We now begin from the two-dimensional Einstein equations and work to construct 
something that resembles the first law. Taking Eq.(3.3) and subtracting half of its 
(two dimensional) trace multiplied by g ab we obtain 

- V 6 V aP + -g ab Up = 47iGp(T ab - T c c g ab ) + 2vGpT c c g ab . (3. 16) 

Contracting both sides with p£ a V b p and using our definition of k and M above we 
find 

-C a V« ^ = CV a M + l -T c c CV a V (3.17) 

where A is the area of the apparent horizon and V = 47rp 3 /3 is the volume it encloses. 
This equation has been noted by Hayward [11] as first law of black-hole dynamics and 
we see here that it also applies to quasi de-Sitter space-times. To zeroth order in e, 
the trace of the stress energy tensor is T 6 C = 2P, where P is the pressure. Therefore, 
it looks natural to identify the last term in the last equation as a pressure term, 
P ~ T c c /2. 

In summary, we identify ^- as the temperature associated with the apparent 
horizon; the mass M and volume V correspond to the mass an volume within the 
horizon. We assume that even if we do not have well-defined expressions for the total 
energy or volume outside the horizon, we still have 

5E out = - 5E m (3.18) 

For 5E in in particular, one can check from the signs in Eq.(3.5) and Eq.(3.7) that 
both the current and the directed surface element are pointing inwards from the 
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horizon. Then, consistently writing all quantities relevant to the space outside the 
horizon (analogous to the expressions for the black hole, where the quantities are 
inside the horizon) Eq.(3.17) can suggestively be written as 

(tf) T- = -CVaM -PC a V a V . (3.19) 

v -yT- ' +SE mt +P5V out 

I ob 

Here to make the analogy more consistent we should have multiplied the equation 
by the normalization factor, B~ l d\, but since in the end this factor appears equally 
in all terms we can ignore it. 

For the unperturbed metric we evaluate this expression for a scalar field source, 
taking 5E in = —5E out to be defined from Eq.(3.5). We find 

2M*H 2 e = $ . (3.20) 

Interestingly, although the potential V(0o) appears in each of the terms on the right 
hand side of the thermodynamic equation, all dependence on V(4>q) cancels out in 
arriving at the expression above. The usual derivation of Eq.(3.20) in GR uses the 
continuity equation combined with the (00) Einstein equation (the first Friedmann 
equation). The thermodynamic version basically inverts that reasoning: in addition 
to the equation above, if we assume we still have the continuity equation T a g. a = 0, 
then the (3 — equation relates the time derivative of energy density pe to energy 
density and pressure. Combining that with Eq.(3.20) above, we recover the time 
derivative of the (00) Einstein equation: 

p E = -3H(p E + P) (3.21) 



, 2 r + V(<Po)J = SH^ = 6MpHH 

So we "recover" the Einstein equation up to an integration constant, similarly to 
Jacobson's result [6]. 

The result in Eq.(3.20) is the same as found by Frolov and Kofman. To match 
their starting expression, however, consider the vector £ a defined in Eq.(3.4), the 
radial generator of p = constant surfaces. For exact de Sitter, £ a becomes null and 
parallel to ( a . Using Eq.(3.11) it can be shown that, on the horizon, C a V a M = 
—i.Tip^ a T ab V b p. Then we can write Eq.(3.17) as 

CV a (^j E£ = -A7rpCT ab V b p + \TX a V a V . (3.22) 

Finally, substituting for T c c using the trace of Eq.(3.3) we find the simple form pre- 
viously derived in [16], 

-CV Q A = 16nGAT b a CV b p (3.23) 

which interestingly relates the change of area of the apparent horizon with the flux 
of energy passing through the p = constant surface. 
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4. Including Perturbations 

To examine what happens in the context of standard slow-roll inflation or eternal 
inflation of the stochastic type we need to include fluctuations in the analysis above. 
Since we have so far used only the Einstein equations and geometric identities, we ex- 
pect the "first law" in Eq.(3.19) to hold whenever the Einstein equations do. We will 
explicitly demonstrate that here, including fluctuations of type required for stochastic 
eternal inflation. 

4.1 Perturbing the spherically symmetric metric 

Since we have used spherical symmetry to establish a framework for our calculations, 
it is convenient to discuss metric perturbations in a format that shows explicitly 
when that symmetry is preserved. The full theory of metric perturbations about a 
spherically symmetric background is well-studied because of its application to black 
holes, and we will follow that literature (concisely reviewed for example in [35]). 

All metric fluctuations, which need not themselves preserve the spherical sym- 
metry, can be written as expansions in spherical harmonics multiplied by functions of 
r and t. However, to illustrate our main point we don't need to consider all possible 
modes. To see why write the the source term, fluctuations in the scalar field, as an 
expansion in spherical harmonics [16]: 



where je(kr) are spherical Bessel functions, Y^ m (#, (j>) are spherical harmonics, H^ l \kr) 
are Hankel functions of the first kind, r ~ —l/(aH) is conformal time and the cre- 
ation and annihilation operators aum, %im ac ^ on the Bunch-Davies vacuum. In 
other words, this is the usual solution but with momentum expanded in spherical 
harmonics rather than Cartesian modes. In particular, the modes still have the be- 
havior of oscillating while well inside the horizon (kr 3> 1) and freezing out on scales 
of the horizon size. For i ^ modes, je{kr) — > as kr = kr — >■ and even on 
horizon scales, r ~ (aH)~ l , these modes are less important. We need only consider 
the £ = mode which dominates there. 

This makes life simpler in terms of the metric perturbations, as we can consider 
only the I = mode perturbations (which preserve the spherical symmetry). For a 
generic spherically symmetric metric, these are [35] 




(4.1) 



ds 2 = ( 9ij + hijY 00 ) dx i dx j 

+p 2 {x i ) laP {x a ){\ + Y 00 K)dx a dxP . 



(4.2) 
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We can write this in a form that is more familiar by defining 

-2$ = h tt Y 00 (4.3) 
-aB r = h rt Y 00 
-2a 2 ^! = h rr Y 00 
-2a 2 r 2 V 2 = r 2 Y 00 K 

and applying this to the usual FRW metric: 

,2 / 1 i O^F>\^2 Q„ D I „2/l OiTr N J„2 i „2„2/l OiTr \f^o2 i „;„2/lJi2\ 



= -(1 + 2<$>)dt 2 - 2aB r drdt + a\l - 2V 1 )dr 2 + aV(l - 2^ 2 )(d8 2 + sin 2 

(4.4) 

A generic gauge transformation on this mode takes x a — > x a = x a + E a with E a = 
y 00 (^o, £i, 0, 0). Our perturbations transform as 

$ -> $ + i , fl r fl r - a- 1 ^ - d r £o) + 2-6 (4.5) 

a 

-f + a~ 2 <9 r £i , ^2 -> ^2 + - "ft 

a err a 

For this limited case, we cannot define gauge-invariant variables, but we can choose 
a gauge that looks familiar by setting B r = and = ^2- Then comparing the 
Gix an d G22 Einstein equations we find $ = ^ and so we have the familiar conformal 
Newtonian gauge (also used by Frolov and Kofman): 

ds 2 = -(1 + 2<&)dt 2 + (1 - 2$)a 2 (dr 2 + r 2 dtt 2 ) (4.6) 

For the perturbed metric the apparent horizon is at 

p H = H-\l + $ + H- 1 ® - (aH)- 1 ®^) . (4.7) 

The (0, 1) Einstein equation gives 

($ + H$) }r = -47rG(-0 o ^),r (4.8) 

and the (0, 0) Einstein equation gives (using the previous equation and the equation 
of motion) 

1 ) $ = ^ ( ^ ) . (4.9) 



V47rG'a 2 2 / dt\d>J 



Substituting these expressions into 2<fi5(f) = (p 2 (5<fi /</>)" + {<j>S(j))' and adding the top line 
of Eq.(4.8) gives an expression we will be able to compare to the 'thermodynamic' 
equation: 

$--$ r +4r$rr + #$+4 ( --a)$ r = 87TGty[ 6<f>--6<f> r ) . (4.10) 
a ' a 2 ' a 2 \r I ' \ a 
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4.2 Comparing the 'thermodynamic' equation 

Now we will compare the result of the Einstein equations above to the 'thermo- 
dynamic' equation. We again evaluate Eq.(3.19), including all terms first order in 
metric perturbations and in scalar field fluctuations. We find (for arbitrary values of 

$ - 2a~ 1 $ r + cT 2 $ rT + = 8irG<j> (5(f> - a _1 50, r ) • (4.11) 

This agrees with the Einstein equation above, Eq.(4.10), up to a term which vanishes 
on the horizon. (It also agrees with the final expression found by Frolov and Kofman). 
Near the original horizon, rH Q = 1 and the expressions are equal. The position of 
the horizon has been perturbed by an amount $ so that $ < 1 is required. Using the 
second line of Eq.(4.8), the condition $ < 1 implies H^fe < M p , which is obviously 
satisfied during inflation. 



5. Eternal Inflation 

Despite the serious difficulties presented by the global picture of an eternally inflating 
multi-verse, the local picture of a Hubble patch in which the inflaton field fluctuates 
up the potential has a well-controlled perturbative description in the sense that 
metric and field fluctuations are small. In some gauges, for example a slicing where 
we choose to follow constant hypersurfaces, the fluctuations may not look small, 
but this is no obstacle to choosing a different slicing in which to follow the story 
[36, 37]. 

To illustrate this, we first verify that in the gauge above the metric fluctuations 
are small even in the eternally inflating regime. The condition for eternal inflation 
is that quantum fluctuations, which may move the inflaton field up or down along 
the potential, are as large on average as the classical distance the field would move 
down the potential in a Hubble time. In other words: 



m ~ i (5.D 



^0 

H 

2e ~ — 
M, 



V 



Since the curvature perturbation which remains constant outside the horizon has an 
amplitude 

V / (m = T7^ (5-2) 



M p V2e 



we see that at the eternal inflation boundary on slices of constant density (5<f) = 0), 
curvature fluctuations are of order one. However, this does not cause any problem 
with our analysis in the previous section. There, we saw that our solution for $ 
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was small regardless of the condition in Eq.(5.1). In other words, we found that $ 
evaluated using the horizon crossing expectation value for S(j) was 

m^E^l.J^ (5.3) 
' ' |0o| M*Hi [ ] 

so that even if the first factor is order one, the second factor must be small while we 
have inflation (and is smaller the flatter the potential). 

Still, something interesting is happening in the eternal regime. Thanks to the 
quantum fluctuations H can remain nearly constant as required but its derivative, 
H, is no longer determined by the classical motion of the scalar field. In the usual 
GR calculation, we take a linear combination of the (00) and (11) Einstein equations 
to solve for H, and it is exactly this linear combination that the 'thermodynamic' 
version of the equation cares about (we've asked for how the horizon area, essentially 
determined by H, is changing). So, in fact the thing that is breaking down in the 
eternally inflating regime is the separation into background and fluctuation we did 
in the previous sections, writing separately Eq.(3.20) and Eq.(4.11). Using the naive 
scaling that in the eternally inflating regime 5<p ~ H8(j) ~ </>o, we see by comparing 
the right hand sides of Eq.(3.20) and Eq.(4.11) that the source terms for H and $ are 
roughly the same size (and there are other terms like (<5</>) 2 , which we dropped above, 
that are also of the same size). This problem is the same in both the standard GR 
treatment of the problem and in the thermodynamic treatment since we are merely 
looking at a particular combination of the Einstein equations restricted to a particular 
surface. 

In practice, we don't know how to formally treat a series of such steps, and so the 
stochastic inflation picture [38] was developed to capture the long-time behavior of 
the system. This has provided a framework for thinking about a more global picture 
of eternal inflation, where we ask about statistics of the field on very large scales, 
coarse-grained over roughly Hubble scales. However, if we wanted to go further with 
the local picture, the calculation is in some ways similar to what must be done for 
an evaporating black hole: we want to calculate the back-reaction of a quantum 
process on the metric. Some difference (and difficulty) comes from the fact that 
we must consider a system that already has classical evolution (e ^ 0) to see the 
quantum effects for de Sitter, and so the back-reaction calculation is in practice 
more complex. (On the other hand, there has been some previous work applying 
the stochastic formulation to black hole evaporation [39]). However, as with the 
evaporating black hole, we may not need to actually do that calculation precisely in 
order to learn something from the system. 

Of course there is also a question of the second law, and whether perhaps there 
is a generalized second law for de Sitter space by analogy with the evaporating black 
hole. A traditional approach to this problem has been to rely on the fact that for any 
given Hubble patch the decrease in de Sitter entropy due to stochastically increasing 
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H is small [40] since the patch eventually exits eternal inflation back into the slow- 
roll regime. It would be very interesting to investigate if something further can be 
said. 

6. Discussion 

In summary, we have revisited the question of thermodynamics and eternal inflation 
of the stochastic type. We have found that a version of the first law, suggested 
previously in [11], but largely examined in the context of black holes, seems to be 
reasonably well- motivated and consistent for quasi-de Sitter space as well. Using 
this expression, we demonstrated that a local description of eternal inflation of the 
stochastic type falls within the regime where the thermodynamic picture can be de- 
fined. This is not surprising, but is a consequence of the fact that there is a pertur- 
bative, local description of stochastic eternal inflation in GR and we have used only 
the Einstein equations and geometric identities to arrive at the thermodynamic equa- 
tion. A more non-trivial test of the suggested "first law" would come from checking 
the next order corrections in e, which would help determine the best definitions for 
surface gravity and entropy. It may also be useful to consider if the thermodynamic 
relations change when one allows perturbative departures from the Einstein Hilbert 
action (higher curvature terms). For black holes, this has been studied since the 
1980's. In the context of Jacobson's work 1 , the investigation began in [41] and more 
recent work includes [42, 43, 44, 45]. Related considerations for static, spherically 
symmetric spacetimes are considered in [46] and FRW spacetimes in [47]. 

However, the implications of even this lowest order result are interesting: both 
black hole evaporation and stochastic eternal inflation are clues to infra-red puzzles 
in gravity. There are hints that the same long time scale that sets the black hole 
evaporation time indicates where an effective field theory description first breaks 
down [5], t ev = RS. Something similar appears in de Sitter computations (for a 
recent summary see [48] and references therein) where R and S are the radius and 
entropy of the (quasi) de Sitter space rather than of the black hole. Those hints 
suggest that, in the regime of dynamical gravity, the global de Sitter solution is just 
an approximation; more generally, it seems that the validity of semi-classical GR 
may be limited in the infra-red 2 , and not only in the UV as naively expected. One 
of the more interesting claims of the thermodynamic approach is that it gives insight 
into gravity independent of the particular microphysical details of the theory. As 
discussed in the introduction in some detail, we expect the regime of validity of the 

1 We thank the referee for bringing this literature to our attention. 

2 Othcr research directions, somewhat orthogonal to this discussion, share similar conclusions. 
An infra-red geometric modification of semiclassical GR has been envisioned [49, 50] such that 
(Tap) evaluates the same as in flat space and the non-local terms responsible for the backreaction 
are absent from the beginning. 
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thermodynamic treatment to be the same as that of semi-classical gravity. If this 
idea has depth, stochastic eternal inflation is an ideal laboratory, complementary to 
the evaporating black hole, in which to test it. 

Alternatively, the thermodynamic analogies may be indications that we should 
not try to quantize gravity in the traditional way. This point of view would be 
strengthened if some thermodynamic argument could provide testable effects beyond 
Einstein gravity. It is tempting to speculate further and give an example inspired 
by the present calculation. As in all the thermodynamic versions of the Einstein 
equations derived so far (e.g. [6, 16, 8, 19]), our first law (3.19) is oblivious to a 
cosmological constant term in the energy momentum tensor. This can be checked 
explicitly by substituting T a f, = —Kg a b on the RHS of the equation and noting that, in 
this case, P = —A. That is, in a cosmological setup, our equation (together with the 
continuity equation) is equivalent to the time derivative of the Friedman equation, 
and A is just an integration constant. At present, however, this is no more than a 
curiosity that appears thanks to our choice of definition of SE. If gravity were really 
insensitive to A, this would suggest deviations from our current picture of GR, such 
that the quantum cosmological constant generated from matter fields is zero, or that 
the value is dynamically irrelevant. For example there should exist a mechanism 
that, during a phase transition, adiabatically readjusts to zero the "zero mode" 
vacuum energy contribution from the matter fields. Because they "just" reproduce 
the Einstein equations, thermodynamic analogies have not revealed a mechanism of 
this type so far, but this could be yet another direction to explore. 
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A. Surface gravity 

A.l Review of surface gravity for a static spherically symmetric space- 
time 

Consider a static spherically symmetric space-time with metric 3 

3 This is not the most general static spherically symmetric metric but for perfect fluid sources 
this is sufficient. For more general cases, see the literature on "dirty black hole" metric solutions. 
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ds 2 = -f{ P )dt 2 + fipy'dp 2 + P 2 dn 2 . (a.i) 

We now follow [51], keeping track of the signs carefully. If a particle is held stationary 
at constant p its acceleration, A a can be computed as 

A a = u°pu , (A.2) 

where u a = dx a jdr is the four velocity along this trajectory and r is proper time. 
One can easily see that A = 0, while A p = |/'(p). Therefore, the magnitude of the 
local force required to hold a unit mass at this position (equated to the magnitude 
of acceleration) is 

A= l -r ll2 { P )\f\ P )\. (A.3) 

Note that at the horizon f{pH) — and the force is divergent. However, if this mass 
is held at that point with a string by an observer at a different radius p, energy 
conservation arguments [52, 51] show that the force that she exerts on the string is 
in fact multiplied by a redshift factor, 

A=(M) 1/2 A. (AA) 



J(P), 

This force is generally finite as the infinite force at horizon is diluted by an infinite 
redshift. By choosing f(p) = 1 we arrive at the expression for the surface gravity 

*=\\f{pH)\. (A.5) 

Note that the above definition is intrinsically non local as it relies on the position 
of the reference observer at p for which f(p) = 1. The obvious choices for a black 
hole and de Sitter in static coordinates are to put the observer at p — > oo and p = 
respectively. 

One can also make use of time translation Killing vector to define surface gravity. 
For a static metric £ a = dx a /dt is the time translation Killing vector. Note that 
any constant multiple of £ a is also a Killing vector. However, £ Q is conventionally 
taken such that £ a = f l ^ 2 u a and at f(p) = 1 it normalizes to £, a £, a = —1. On the 
other hand, at the horizon £ Q is null and both tangential and normal to the three 
dimensional horizon itself. Thus, one can also identify this hyper-surface as the place 
where = 0. Since £ Q remains tangent to the null generator of the horizon it also 
satisfies the geodesic equation for some constant k: 

U-A* = (A.6) 

which matches the k introduced in Eq.(A.5) as we now review. First making use of 
the Killing equation, 

Up + = 0, (A.7) 
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equation (A. 6) can be turned into: 



"foa^ = ( A - 8 ) 

or we can write it as 

-(£%,);« = 2k& (A.9) 

which is not surprising since at the horizon £ a is normal to the — (C^Cp) constant sur- 
face. Now to calculate different sides of this equation and obtain k we use Eddington- 
Finkelstein ingoing or outgoing coordinates depending on which one is well suited 
to describe the horizon. For example, if the outgoing null ray is tangential to the 
horizon (as for the Schwarzschild metric) then the ingoing coordinate is appropriate: 

ds 2 = -f(p)dv 2 + 2dvdp + p 2 dtt 2 . (A. 10) 

However, if the ingoing null ray is tangential to the horizon (as for the static de-Sitter 
metric) then the outgoing coordinate is appropriate: 

ds 2 = -f(p)du 2 - 2dudp + p 2 dtt 2 . (A. 11) 

For Schwarzschild, u constant and p constant surfaces coincide at black hole horizon 
and are degenerate while for de Sitter v constant and p constant surfaces coincide. 
Transferring £ Q to this coordinate and calculating £ Q for these cases yields £o — — f 
and £i = ±1 where the plus and negative signs correspond to v and u coordinates 
respectively. This implies that 

= f,c 

= f'p,a (A. 12) 

and also that at the limit where / = 0, 

£ Q = ±p, Q (at p = p H ) (A. 13) 

and together with Eq.(A.9) we finally obtain: 

« = ±\f'{pH). (A.14) 

It is easy to see intuitively why the signs here agree with \f'(pH)\ i n Eq.(A.5). The 
positive (negative) sign is for the case of outgoing (ingoing) null ray being tangential 
to the horizon, so it must be blocking the region p < Ph{p > Ph) from observers. 
Therefore the observer p must be at p > pn (p < Ph)- Assuming f(p) is monotonic 
and given that }\ph) = and f(p) = 1, we expect f\pn) > (/'(ph) < 0). 
Notice that it is now evident why for these metrics defining horizon as (Vp) 2 = is 
consistent with = 0. 
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A. 2 Surface gravity in time dependent spherically symmetric spacetimes 

Let us start by writing a general spherically symmetric metric in the form 

ds 2 = g ab (r,t)dx a dx b + p 2 (r,t)dtt 2 . (A. 15) 

Here g ab is the metric for the (t, r) coordinates and is sometimes referred to as the 
base metric and p 2 is the induced metric of two spheres at constant r and t and 
related to their areas through p 2 = A/Att. A marginally trapped sphere or apparent 
horizon is defined such that at least one of 9 + = d\nA + /d\ or = d\nA_/d\ 
vanishes along the null ray itself and consequently (as was explained in section 2) 
(Vp) 2 = at apparent horizon. If 9 + and are zero at the same time, this is a 
bifurcation surface, but for the interest of this paper we consider the case in which 
only one of them is zero and it changes sign from one side to the other side. The union 
of all marginally trapped surfaces formed through extending the apparent horizon 
into past and future is called a trapping horizon and we will refer to it as E. Spheres 
for which + 0_ < 0, are called untrapped surfaces and one can fix the orientation of 
ingoing and outgoig in the region they reside by taking 6L < and 9 + > [21]. On 
the other hand the region on the other side of E is where trapped surfaces, defined 
by 9 + 9_ > 0, reside. To see how trapping horizons in cosmology differ from trapping 
horizons for blackhole and whitehole see table A. 2. Next we consider cases where 
the apparent horizon is either a future outer trapping horizon or past inner trapping 
horizon and by that we mean: if = (6L = 0) then 9^d~9 + > (6* + <9 + 6L > 0). In 
the static limit this includes cases where apparent horizon is protecting the trapped 
side from sending light rays to untrapped side. 





Blackhole 


Whitehole 


Expanding FRW 


Contracting FRW 


Untrapped 


9+ > 


9+ > 


9+ > 


9+ > 


spheres 


9_ < 


9_ < 


9_ < 


9_ < 




(outside) 


(outside) 


(inside) 


(inside) 


Trapped 


9+ < 


9+ > 


9+ > 


9+ < 


spheres 


9_ < 


9_ > 


9_ > 


9_ < 




(inside) 


(inside) 


(outside) 


(outside) 


Marginally 


9+ = 


9_ = 


9_ = 


9+ = 


trapped 


<9_ < O(future) 


9 + > 0(past) 


6+ > 0(past) 


6L < O(future) 


sphere 


<9_6> + < (outer) 


d + 9- < (outer) 


d + 9_ > O(inner) 


d-9 + > 0( inner) 



Table 1: Characterizing trapping horizons based on expansion of congruence of light rays. 
The first row fixes the orientation of ingoing and outgoing vectors in untrapped region and 
removes the ambiguity for trapped region. d± denotes coordinate derivative along l±. 

As was discussed earlier, the definitions of k are intrinsically nondocal and they 
also rely on the fact that the horizon is null and £ Q is a Killing vector that satisfies 
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the Killing equation (A. 7). These criteria lead to 



£ a (£/J;a + =0 



(A.16) 
(A.17) 



and hence any linear combination of above equations can be used to define surface 
gravity as long as normalization of £ Q is fixed. Note that also £ a becomes null at 
horizon and it is tangential to the null ray whose expansion vanishes. The Kodama 
vector £ a = —e ab 'S7bP is also null at marginally trapped surfaces and it is tangential to 
the null ray whose expansion vanishes (which we review in a few lines). Moreover, it 
goes into the Killing vector in the stationary limit of metric type Eq.(A.l). It seems 
natural then to ask if substituting Kodama vector for Killing vector in the definition 
of surface gravity is possible. 

Unfortunately since £ Q is no longer necessarily a Killing vector, the LHS of equa- 
tion Eq.( A. 17) remains proportional to £p but the RHS of Eq.(A.16) is not necessarily 
zero and in general it is not even proportional to Even when £ is a null hyper- 
surface and the RHS of equation Eq.(A.16) is proportional to different linear 
combinations of the two equations result in different definitions for k. However, 
as Hayward suggested [11] since Eq.(A.17) remains valid for spherically symmetric 
metrics one can use that for defining surface gravity. Of course now the question of 
normalization remains since k clearly changes if £ is rescaled by a constant. How can 
we determine that the normalization of £ is consistent with normalization of Killing 
vector for null horizons? We review some proposed solutions in what follows. First, 
we show Kodama vector is in fact tangential to a appropriate null ray at S, second 
that its normalization is consistent with the local normalization of Killing vectors 
and last that this definition is consistent with 



We denote the outgoing and ingoing null vector as I" and I" respectively. Solving 
the second order null equation for this region yields: 



To determine which solution corresponds to which null vector, r constant surfaces 
were taken to be time-like (except at horizon itself) and t constant surfaces to be 
space-like or null, thus <?oo must not be positive and gu must not be negative. On 
the other hand the condition 9 + = 0(6L = 0) is equivalent to Z+V a p = 0(Z" V a p = 0). 
Therefore, 




(A.18) 




(A.19) 



V a p OC l± 



a 



for 6± = 



(A.20) 
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In fact, substituting Eq.(A.19) in Eq.(A.20), the indices of Kodama vector can be 
lowered in a consistent way: 



£ a = ±V a p for 9 ± = (A.21) 
oc l ±a for 9± = (A. 22) 

As we discussed before, even though Eq.(A.6) is a local equation that can be applied 
to any Killing horizon it does not provide a local definition of surface gravity, since £ 
was not normalized locally, independent of p. To overcome this challenge an alterna- 
tive method has been put forward [22] which relies on first normalizing the auxiliary 
null vector field N a , defined by N a £ a — — 1. In fact this method can also be used to 
fix the rescaling of the Kodama vector. The prescription is that one can impose the 
normalization by fixing the nonzero expansion of null rays at X itself: 

9 T = T~ p for 8± = (A.23) 



For example let us rescale £ Q by a constant £ a — > £ a = ±C V a p. Associating N a to 

dx a I 



the other null beam tangent to the null geodesic, N a = ^tt-\n , the above condition 



leads to: 

e T = - P , a N a 
p 

= T~ p ■ (A.24) 

Thus p, a N a = =Fl and further fixing of N a ^ a = — 1 immediately results in C = 1. 
Now that £ is normalized it is justified to use Equation (A. 17) to calculate k. This 
equation can be expressed as 

£ a V [a £ 6] = k£ 6 . (A.25) 

To prove this equation note that V[ a £b] = [afo]V[o£i], where [ab] represents the per- 
mutation symbol. This implies 

£ a V [a £ b] = -HV [0 £i]e ac V c p 
_ V [0 £i] , 



-Vbp 
-9 

V [0 £ 



m , (A.26) 



V~9 

where to arrive at the last line we used Eq.(A.21). We now see that Equations (A. 17) 
and (A.25) hold at S if 



Substituting for £ a = g ac ^ c and using the definition of Kodama vector, we observe 
that 




ge cd V c (g de e e N f p) 



(A.28) 
(A.29) 
(A.30) 



where in the last line we have used e cd e e fgd e = g and that Levi-Civiata tensor is 
covariantly constant with respect to (1 + 1) covariant derivative. Therefore 



It is worth noting that, for a slow-roll inflationary back-ground, the definition put for- 
ward in [33] (based on projected acceleration of the vector tangential to the trapped 
surface) differs from the above at 0(e). 
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